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Q^' We make quantitative predictions for the rapidity and centrality dependencies of 

liadron multiplicities in AA, pA and pp collisions at the LHC energies basing on the 

• • . ideas of parton saturation in the Color Glass Condensate. 
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X 

c5 ■ INTRODUCTION 

At high energies QCD is expected to enter the new phase : the Color Glass Condensate 

eading to parton saturation 



S]GC) which is characterized by strong coherent gluon fields leading to parto: 
; 0; 0; Ol • Prcvlously, we have applied this approach 11 , [l^ 13 , 14 , [3| 



to describe 

the wealth of experimental data |6|, 13, |9|, llO| from RHIC. The LHC will allow to extend 
further the investigations of QCD in the regime of high parton density. This is because 
the new scale of the problem, the saturation momentum Qg, will become so large {Ql ~ 
5 — 10 GeV^) that a separation of CGC physics from non-perturbative effects should become 
easier. The main objective of this paper is to give predictions for the global characteristics of 
the inelastic events in nucleus-nucleus, proton-nucleus and proton-proton collisions at LHC 
energies basing on the ideas of parton saturation in the Color Glass Condensate (CGC). 
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To understand better the differences implied by a higher energy of the LHC, let us start 
with the main assumptions of the approach we used to describe the data from RHIC: 



1. At Bjorken x < 10~^ the inclusive production of partons (gluons and quarks) is 
driven by parton saturation in strong gluon fields as given by McLerran-Venugopalan 
model 



by ] 

1 



2. The region of x ~ 10~^ (accessible at forward rapidities at RHIC) is considered 
as the low x region in which as ln(l/x) ^ 1 so the quantum evolution becomes 
important; we assume that ^ 1 to keep the calculation simple and transparent; 

3. We assume that the interaction in the final state does not change significantly the 
multiplicities of partons resulting from the early stages of the process; this may be 
a consequence of local parton hadron duality, or of the entropy conservation. There- 
fore multiplicity measurements are extremely important for uncovering the reaction 
dynamics. However, we would like to state clearly that we do not claim that the 
interactions in the final state are unimportant. Rather, we consider the CGC as the 
initial condition for the subsequent evolution of the system, which can be described 
: or example by means of hydrodynamics (such an approach has been followed in Refs. 

Even a superficial glance at these three assumptions reveals that the conditions for the 
applicability of our approach at the LHC improve. Indeed, at LHC energies the value of x will 



be two orders of magnitude 



methods of low x physics 



ower than at RHIC. This makes the use of the well-developed 



12 



3,0,0 



22] better justified. At LHC energies we have 



a theoretical tool to deal with the high parton density QCD in the mean field approach ( 
so called Balitsky-Kovchegov non-linear equation 4]), or on a general basis of the JIMWLK 
equation j^; even more gene ral approaches may be possible (see for example the lancu- 



Mueller factorization 



2G 



2l|). However, despite a number of well developed approaches 
which could be applied at low x we would like to warn that even the LHC energy is not 
high enough to apply any of the methods mentioned above without discussing possible "pre- 
asymptotic" corrections to them. 

Consider for example the determination of the value of the saturation momentum - the 



key scale in the CGC phase of QCD. As was noticed first in Ref. 



2J] the value of the 
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saturation scale is affected by the next-to-leading order corrections to the BFKL kernel 
which were neglected in all of the discussed above approaches. Their numerical significance 
is so large that they cannot be neglected: if the next-to-leading order BFKL kernel is used, 
because of a large energy extrapolation interval to the LHC the value of turns out to be 
5-10 times smaller than if one uses the leading order kernel ( see detailed discussion in Ref. 
[2^). However the good news is that the NLO corrections appear under theoretical control 
and we can take them into account. 

The paper is organized as follows. In the second section we discuss the geometry of 
nucleus-nucleus and hadron-nucleus collisions and introduce the Glauber formalism we use. 
In the third section, we review the general formalism which we use to evaluate the mul- 
tiplicities; we also discuss the influence of higher order corrections and the effects of the 
running coupling constant on the results. In the fourth section we list the parameters of 
our approach and justify the values we use; we then give a complete set of predictions for 
hadron multiplicities at the LHC energies in Pb — Pb, p — Pb, and pp collisions, including 
the dependences on rapidity and centrality. We then summarize our results. 



THE GEOMETRY OF NUCLEUS-NUCLEUS AND HADRON-NUCLEUS 
COLLISIONS AND THE GLAUBER APPROACH 

At high energies the paths of the colliding nucleons can be approximated by straight 
lines, since in a typical interaction t/s <^ 1 and the typical scattering angle is small. This is 
the most important approximation underlying the Glauber approach to nuclear interactions. 
Other approximations which simplify calculations but are in principle unnecessary are the 
smallness of the nucleon-nucleon interaction radius compared to the typical nuclear size, and 
the neglect of the real part of the A^A^ scattering amplitude. Many quantities characterizing 
the geometry of the collision can be readily computed in this approach; a complete set of 
the relevant formulae can be found e.g. in |27] and we will not reproduce all of them here. 

It is customary and convenient to parameterize the centrality of the collision in terms 
of the "number of participants" Npart ^ the number of nucleons which underwent at least 
one inelastic collision. This number can be directly measured experimentally (at least in 
principle) by detecting in the forward rapidity region the number of "spectator" nucleons 
Ngpect which did not take part in any inelastic collisions; obviously, for a nucleus with mass 
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number A, Npart = A- N^pect- 

The number of participating nucleons in a nucleus- A-nucleus-B interaction de pen ds on 
the impact parameter h. In the eikonal approximation it can be evaluated as (see (2fi|): 

+ B j d\sTB{h-s){l-[l-a,nTA{s)]^] , (1) 

with the usual definition for the nuclear thickness function Ta{s) = dzpA{z, s), normal- 
ized as / (PsTa{s) = 1; is the proton-proton inelastic cross-section without diffractive 
component. For the LHC energies we assumed = 70 mb ([28]). 

From Eq. ((H) the definition of the local density of participants np^^{h,s) is evident; we 
will define its average over the transverse plane as 



In the following we will need to use the average number of participants computed separately 
for nucleus- A and nucleus-B; it is given by 



AB \ /LN _ ^ "-part.Al^, s) ^part(b, s) 



/n^B WTX _ J " part,A\ ^ i party ^ i /o\ 

Obviously, one has for their sum 

where ^part^Al^^ and ^part.Bl^) the integrands of the first term and second term in 
the r.h.s of Eq. ((T)) respectively. 

In table El we give the number of participants and their density (respectively Eqs. ((H) and 
((21)) for Pb-Pb collisions at LHC. 



The corresponding formulae for the proton-nucleus j)A interaction can be deduced by 
setting B = 1 and using a delta-function for the proton thickness function (in the point-like 
approximation for the size of the proton) . We get from Eq. (^ : 

N^artib) = Aa^nTAib) + {l " [1 " a,„T^(6)]^} = Aa,„T^(6) + {l - P^^ib)} . (4) 
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6 

(fm) 


part 


part 

(fm-2) 


b 

(fm) 


part 


r,AB 
' ''part 

(fm-2) 


0.00 


406.9 


2.98 


8.00 


166.8 


2.21 


1.00 


402.4 


2.97 


9.00 


127.5 


1.97 


2.00 


387.8 


2.93 


10.00 


91.9 


1.69 


3.00 


363.2 


2.88 


11.00 


61.1 


1.35 


4.00 


330.3 


2.80 


12.00 


36.2 


0.98 


5.00 


291.9 


2.70 


13.00 


18.3 


0.59 


6.00 


250.6 


2.57 


14.00 


7.5 


0.27 


7.00 


208.3 


2.41 


15.00 


2.5 


0.09 



TABLE I: Mean number of participants and their average density in Pb-Pb collisions at LHC as a 
function of b 

In the previous formula the function PQ^^^b) is the probability of no interaction in a p-A 
collision at impact parameter b; the integration of [1 — PQ^{b)] over b gives the inelastic 
proton-nucleus cross section apA- 

The average number of participants in a p-A collision can be obtained as: 



Jd'b[l-Po{b)] 



(5) 



the first term in the r.h.s. gives the mean number of participants (iVp^^^^) in the nucleus. As 
in the case of nucleus-nucleus collision, we will need to compute the density of participants 
in nucleus A, defined as: 

\ _ {^part,A} _ A 



{n: 



part, A/ 



O^in (^pA 

In practice, the information about the impact parameter dependence is extracted by 
analyzing the data in various centrality bins. The physical observable most frequently used 
to estimate the centrahty of the coUision is the multiphcity of charged particles Nch- We 
will assume that the average value of Nch produced in a collision at impact parameter b is 
determined by the number of participating nucleons N^artip). The actual multiplicity will 
fluctuate around its mean value according to: 



V{Nch,{Nch{b))) 



^27ra{N,„{b)) 



C{{Nch{b)))exp 



MNch{h)) 



(7) 
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where the factor C{N) = 2/[l + erf (^N/ 2a)] is introduced to ensure that the fluctuation 
function V{Nch,N) satisfies Jl^ dNchVlNch, N) = 1. The numerical value of C{N) is 1 
with very good accuracy for almost all cases of practical interest (it can exceed 1 for very 
peripheral collisions, where the number of participants and consequently Nch is small: in 
such a case it is important to include the factor C{N) to have a correct normalization). 

The parameter a gives the width of the fluctuations: its value is dependent on the ex- 
perimental apparatus, therefore it is not possible for us to predict its value for the LHC 
experiments. For the experiments at SPS and RHIC the value of a varies from 0.5 to 1.5-2. 
We will assume a = 0.5 in the following; uncertainty in this parameter can affect the cen- 
trality dependence of our results. In the case of PbPb collisions, we estimate the resulting 
uncertainty in the density of participants (and thus in the saturation scale, see below) to be 
about 5%; in the case of pPb collisions, this uncertainty can reach 10 -r- 15% for peripheral 
collisions. 

We will also assume the proportionality between N^h and Npart when computing the dif- 
ferential inelastic cross section; this proportionality is not exact, but the shape of minimum 
bias distribution of events which is normally used to fix the parameter a (and the propor- 
ionality constant between Nch and Npart) has been found insensitive to this assumption (see 



The minimum bias differential cross section can be obtained as {N{b) = qNpartip), where 
g is a constant): 



here Po{b) is the probability of no interaction at the impact parameter b: for a nucleus- 
nucleus collision Poib) = [1 — ainTAsib)]^^ where Tab is the overlap function : TAsip) = 
J d'^sTA{s)TB(h — s); in the case of B=l, Poib) reduces to P^^^ib) defined above. In the 
following, all of the formulae will refer to A-B collisions; with obvious modifications they 
are valid also in the p-A case. 

The total nucleus-nucleus cross section is then obtained by integrating Eq. (jH)) over dNch- 



The mean value of any physical observable O (given in terms of the impact parameter b) 





(8) 
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centr. cut 




(fm-2) 


0-100 % 


103.2 


1.33 


0-6 % 


369.0 


2.89 


0-10 % 


346.6 


2.83 


0-25 % 


274.2 


2.62 


25-50 % 


103.7 


1.75 


50-75 % 


27.0 


0.76 


75-100 % 


3.9 


0.14 


0-50 % 


186.7 


2.17 


50-100 % 


15.7 


0.45 



TABLE II: Mean number of participants and their density in Pb-Pb collisions at LHC for different 
centrality bins 



can be computed as: 



(O) 



1 



dN^H ^ 0{b) 



(10) 



(Tab J dNch 

To obtain the corresponding average for a given centrality cut we have to limit the 
integrations in the previous formula in the appropriate way, for instance the expression: 

{O) 



Oih) 



No dNch 



(11) 



gives the average value of the observable O in the fraction of the total cross section defined 
by the limit Nq. In this work the previous formula has been used to compute the mean 
density of participating nucleons (Eq. (0)) in different centrality bins, as shown in table ITTl 



Table IIIII gives the results of Eq. (0) for the case of p-Pb collisions at LHC energy. The 
corresponding densities are obtained according to Eq. (jHl). 
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centr. cut 




0-100 % 


7.41 


0-20 % 


13.07 


0-50 % 


11.31 


20-50 % 


10.29 


50-100 % 


3.58 



TABLE III: Mean number of participants in p-Pb collisions at LHC for different centrality bins 

THE GENERAL FORMULAE 



Let us discuss the main features of the approach we use to describe the production 
dynamics. As in our previous papers jl^ . Il^l we use the following formula for the 



inclusive production 




Am-. 



AttNc 1 ,,2 / ,2^ / / 1^2^ 



(12) 



' d^p - 1 pf 

where Xi,2 = {Pt/ y/s) exp{^y) and (pAi,A2{x, f^t) is the unintegrated gluon distribution of a 
nucleus ( for the case of the proton one of ipA should be replace by ipp.) This distribution is 
related to the gluon density by 



xG{x, 



dktif{x,kf). 



(13) 



We can compute the multiplicity distribution by integrating Eq. (fT^ over pt, namely, 

dN 1 f ^2 j^^'^ . 
dy S J d^p^ 



(14) 



5* is either the inelastic cross section for the minimum bias multiplicity, or a fraction of it 
corresponding to a specific centrality cut. 



Saturation scale 



Let us define two saturation scales: one for the nucleus Ai and another for the nucleus 
A2. We will see below that even in the case of Ai = A^ the introduction of two saturation 
scales will be useful. It is convenient to introduce two auxiliary variables, namely 
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Qs,miniy, W) = min (^Qs{Ai; xi = ^ e (5s(^2; ^2 = e^)^ ; 
Qs,max{y,W) = max (Qs{Ai;xi = ^e'''^),Qs{A2]X2 = ^e^)^ . (15) 



To understand the physical meaning of these two scales we start with the explicit formula 
for Qs which was suggested in Ref. j22| for the description of HERA data on deep inelastic 
scattering and was successfully used to describe the data from RHIC 
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12 



jiii 

3: 



Xo 



A 



Qii^) = Qli^] (16) 



with the central value of A = 0.288 [23]; the value of A has an uncertainty of 5 — 10%. 
Substituting xi = (Qs/W) and X2 = (Qs/W) e^, where W is the energy of interaction, 
one can see that the energy and rapidity dependence of the saturation scale can be reduced 
to a simple formula 

Ql{Ay,W) = Ql{A;Wo) {^/)^' = Ql{A;Wo) (^)' e^^ . (17) 

In what it follows we will use the notation A for A = A/(l + |A) = 0.252 hoping that it 
will not lead to misunderstanding. 

Using Eq. (fT7|) one can see that for a production of the gluon mini-jet at rapidities y 
there are two different saturation momenta: Q^(y4;?/, PF) and Q^A; —y,W), even for the 
collision of identical nuclei (see Fig. Q). Fig. Q shows that the density is quite different in 
two nuclei since at y ^ (say y > 0) one of the nuclei probed at relatively large x = Xi > X2 
is a rather dilute parton system while the second nucleus has much higher parton density 
than at y = 0. Therefore, for an A + A collision at y > Qs,min = Qs{A; —y,W) while 
Qs,max = Qs{A:, y, W). In the case of a collision of two different nuclei we need to take into 
account the A— dependent values of Qo{A; W) in Eq. (fTTj) . 

The saturation scale is the main parameter of our approach and we need to understand 
clearly the energy dependence of Qg if we want to make predictions for the LHC energies. 
The first basic result on the behavior of this scale is the power-like energy dependence which 



ows directly fror a Q CD for fixed QCD coupling. As was shown in a number of papers 



25 



3, 



31 



32 



33 



3JI the energy dependence of the saturation scale does not depend on the 



details of the behavior of the parton system in the saturation domain but can be determined 




n 



FIG. 1: The CGC approach for nucleus - nucleus collision with the saturation of parton density. 



just by using the perturbative QCD approach in the BFKL region Indeed, consider the 
dipole-target scattering amphtude in the double Melhn transform representation, namely, 

N{y, r') = [ ^ + Hr^^hcn) Niu, 7) . (18) 

The BFKL equation determines the value of u at which N{ij,'~f) has a pole: 

w = asxh) (19) 

with a specific function x which can be found e.g. in Ref. |i2^]; we denote as = ^c^s/tt. 
To find the energy dependence of the saturation scale we first need to find a critical value 



of 7 = 7cr defined by the equation 



33, 



34| 



1 - 7cr d'j 

The meaning of this equation is the following: in the semi-classical approximation (see Ref. 
js^ and references therein) the scattering amplitude N{y,ln{r'^AQfj^)) has the following 
form: 

N{y, e = Hr^Alco)) = const x exp [u;{y, Ov - (1 - liv, 0)^] ■ (21) 

The boundary of the saturation region is determined by the unique (critical) trajectory for 
the non-linear evolution equation in the {y, ^) plane for which the phase Vphase = ^^{y, 0/(1 ~ 
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liVy 0) the group Vgroup = —duj{y, ^ ) / d'y{y, ^ ) velocities are equal. The physical meaning 
of this trajectory can be illustrated by an analogy in geometrical optics: the boundary which 
it defines is similar to the focal refiecting surface (therefore, one can see that the surface of 
the Color Glass shines!). The equality of phase and group velocities thus gives the equation 
for the saturation scale: 

d\n{Q^,{x)/AlcD) _ _ Xilcr) _ . 

— JhTM — - ""T^ = ^ ^^^^ 

For fixed as Eq. ^ leads to 

Ql{x)=Ql{^y (23) 

with A given by Eq. (j221)- The numerical analysis of the value of A can be found in Ref. 
2^. The main conclusion from this analysis is the fact that the value of A is sensitive to 
higher order correction in as- Therefore in this paper we choose to fix the value of A from 
the phenomenological approach, see Eq. (fTTj): we consider Eq. as a justification for the 
use of such a parameterization. 

Another observation on the equation for the saturation scale Eq. ()22|1 is that the value 
of 7cr is stable with respect to higher order corrections and almost does not depend on the 
value of the QCD coupling (see Ref. j3|). This fact helps us to solve Eq. in the case of 
running as- The running of the coupling constant as leads to an additional dependence on 
Qs in the r.h.s. of Eq. (f^^ : from Eq. (j^^ using the explicit form of the running coupling 
constant we find 

dlniW/Wo) P2l-lcrHQ'sm/AlcD) HQ'siW)/AlcDy ^ ' 

as a result, the dependence on QsiW) has become explicit. Integrating Eq. we obtain 



QiiW) = A'qcd exp(^v^2 5 1n(iy/iyo) + hi\Ql{Wo) / AIcd)) . (25) 

where QliWo) is the saturation scale at the energy Wq which we used as an initial condition 
in integrating Eq. ()24|1 . Here as well as in the rest of the paper A^qj^ is defined by as = 
47r//?2 1ii(Q^/Aq^£,) and in numerical applications we took A^^jj^ = 0.04 GeV^ with /?2 = 
11 — 2/3 Nf where Nf = 3 is the number of fermions (number of colors Nc = 3)- We fix the 
value of S through the empirical value of A as given by Eq. ()23|) and the value of saturation 
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scale for the Au nucleus at fixed energy of W = 130 GeV, y = 0, corresponding to the cut 
of — 6% of most central colhsions, Q^q = 2 GeV^, so that 6 = X ln(Qfo/AQ^^). 

The formula Eq. reproduces all general features expected for the case of running 
QCD coupling; in particular, one can see that the saturation scale (1^ does not depend 
on the mass number of the nucleus in the limit of high energies 0, \^'\ - the parton wave 
functions of different nuclei in this limit become universal. It is easy to generalize Eq. ()25|) 
to y 7^ by replacing ln(Vr/Wo) by ln(iy/iyo) + thus we have the following final formula 
for the case of running as'- 

Qliy, W) = Ajc^exp ^^2X\n{Q%/AlaD)[HW/Wo) + y] + \n\Ql{Wo)/AlcD)) ■ (26) 



Formulae for the multiplicities 



To derive the final expressions for the multiplicity it is convenient to re-write Eq. ()14|) 
using the fact that the main contribution to Eq. (jl4j) is given by two regions of integration 
over kt'. kf <^ pt and \pt — kt\ <^ pt] this leads to 

— - - J dp, \E— - - -j;^^ 



X 



dp^ 
Pt 



^A^Xupf] 



Pt 



dkt ^A2{x2,kt) + ^A2{x2,Pt] 



■pt 



dkl ipA^{xi,kf 



S N^-l JO PI 

where we integrated by parts and used Eq. ()13p . In the KLMN treatment 
we assumed a simplified form of xG, namely, 
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(27) 

Q 



xG{x]pI) 



Spl (1 



X] 



SQlix) (1 



Pt < Qs{x) 



Pt > Qs{x) 



(28) 



I c^sm 

where the normalization coefficient k has been determined from the RHIC data on gold-gold 
collisions. We introduce the factor (1 — xY to describe the fact that the gluon density is 

nn 

small at X ^ 1 as described by the quark counting rules |38L 18 9|. 

We have checked that the simplified form of Eq. ()28j) is adequate for the calculations of 
multiplicity since it is dominated by the low momenta region. At high pt and small x, it 
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was shown 14] that the quantum effects of the anomalous dimension could be extremely 
important. However, at moderate values of x the simple form of Eq. (j^Hj) was used to 
calculate the pt spectra in proton-proton and electron-proton collisions in Ref. \^ and the 
results appear very encouraging. 

Having in mind Eq. (j2Hl), let us divide the pt integration in Eq. (fT^ in three different 
regions: 

!• Pt ^ Qs,min 

In this region both parton densities for Ai and A2 are in the saturation region. This 
region of integration gives 

^ cx —SQl^,^ cx — iVp,,,(A) (29) 
dy as as 

where we have used the fact that the number of participants is proportional to SQ"^, 
where S is the area corresponding to a specific centrality cut. 

2. Qs,max ^ Pt ^ Qs,min 

For these values of pt we have saturation regime for the nucleus A2 for all positive 
rapidities while the nucleus Ai is in the normal DGLAP evolution region. Neglect- 
ing anomalous dimension of the gluon density below Qs,max, we have ipAi{xi,k^) oc 
^ S Qs,min/kt which for y > Vc leads to 



as 



dN 1 2 1 Qs,max 1 i,T / /i \ 1 Qs,max /nn\ 
— OC — i„ln— 2 (X — Npart{Ai) lu — ^ (30) 

dy as Qs,min "5 Qs 

This region of integration will give the largest contribution. 

3. Pt ^ Qs,max 

In this region the parton densities in both nuclei are in the DGLAP evolution region. 

Substituting Eq. into Eq. (jTTj) we obtain the following formula 

^ = Const X S QlrniniW,y) \ x (31) 



1 



Qs,m.n^.y) ^ {HQl^AW,y)/Ql^^SW.y)) + 1} (^i _ ^£=1^ 
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One can see two qualitative properties of Eq. (|HT|). For y > and close to the frag- 
mentation region of the nucleus Ai , Qs,min = <5(^i) and the multiplicity is proportional 
to Npart{Ai), while in the fragmentation region of the nucleus ^2(2/ < 0) Qs,min = <5(^2) 
and dN/dy oc Npart{A2). We thus recover some of the features of the phenomenological 
'wounded nucleon' model 12611 . 



PREDICTIONS 



Choice of the phenomenological parameters 

As discussed above our main phenomenological parameter is the saturation momentum. 
An estimate of the value of the saturation momentum can be found from the following 
condition: the probability of interaction in the target (or "the packing factor" of the partonic 
system) is equal to unity. The packing factor can be written in the following form: 

= {N^ - 1) ttR^ = ^ ^ (^^^ 

where a is the cross section for dipole - target interaction (the size of the dipole is about 
1/Q) and p is the (two-dimensional) transverse density of partons inside the target of size 



e (t^ 
III 



R (see e.g. Refs. |ll|,|l^ for details). 

In the case of the nucleon we do not know the value of R or, in other words, we do not 
know the area which is occupied by the gluons {Sn = ttR"^). However, we have enough 
information to claim that this area is less than the area of the nucleon ( R is less than the 
electromagnetic radius of the proton). To substantiate this claim, let us recall for example 
the constituent quark model in which the gluons are distributed in the area determined by 
the small (relative to the size of the nucleon) size of the constituent quark. Having all these 
uncertainties in mind we use the phenomenological Golec-Biernat and Wuesthoff model 2^ 
to fix the value of the saturation moment in the case of the nucleon target. Namely, the value 
of the saturation moment for proton is equal Qs{P] y = 0,W = 200 GeV) = 0.37 GeV^. In 
Ref. j|41j this value of the proton saturation momentum was used to describe the deuteron- 
gold collisions at RHIC energies. 

In the case of the nuclear target the saturation momentum can be found from the ex- 
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pression for the packing factor 

P.F. = apA = apM = QsK^) ^ ^ (33) 

As we have discussed we do not know the last factor {Sn / Sa) and therefore, Eq. (jHHj) cannot 
help us to determine the value of the saturation momentum for a nucleus. We fixed the value 
of the saturation momentum from the description of the RHIC data on the multiplicity in 
gold-gold colhsions, namelvJ^^j^Gold, y = 0, ly = 130 Ge^) = 2.02 GeV^ for the centrality 
cut - 6% (see Refs. HQ for details). 

As far as energy dependence of the saturation scale is concerned, we used Eq. p6|) and 
Eq. (j^mi with A given by the Golec-Biernat and Wuesthoff model (A = 0.252). However, we 
have to admit that the perturbative QCD estimates described above would lead to a larger 
value of A : A ~ 0.37. Such an uncertainty in the value of A leads to an error of about 
12 — 15% in our prediction for the proton-nucleus and nucleus-nucleus collisions at the LHC 
energies. For the proton-proton interaction it could generate an error as big as about 50%. 

In our main formula given by Eq. (jHT|) we have to fix the normalization factor. As 
discussed in Ref. [12] theoretical estimates lead to a value of Const in Eq. ()31|) which 
appears quite close to the value extracted by comparison with the RHIC data. In this 



paper we use the same normalization factor Const as in [ij], namely Const S Ql ,^^^{W = 
200 GeV,y = 0) = 0.615 Npart- We also need to note that the experimental measurement 
are done at fixed pseudo-rapidity rj, not rapidity y; therefore, as discussed in Ref. jl3| we 
have to use the relation between r] and y, and to multiply Eq. by the Jacobian of this 
transformation h{ri,Qs)] see Ij] for details and explicit expressions. 



Proton-proton collisions present an additional problem caused by the deficiency of the 
geometrical interpretation of pp cross section. As mentioned above (see Eq. ()14|l ) we calcu- 
lated the ratio of the total inclusive cross section to the geometrical area of the interaction. 
This ratio is the measured multiplicity in the case of hadron-nucleus and nucleus-nucleus col- 
lisions. In the case of hadron-hadron interaction the multiplicity is the ratio of the inclusive 
cross section divided by the inelastic cross section. As discussed above, for pp interactions 
we do not know the relation between the interaction area and the value of the inelastic cross 
section. To evaluate multiplicities in proton-proton interactions we do the following: (i) 
fix the ratio SN/cTin at W = 200 GeV using the data for dN/dy{y = 0); and (II) assume 
Sn/c^'iti oc I /(Tin far as the energy dependence is concerned. In other words we assume 



16 



that the area does not depend on energy. The energy dependence of the inelastic cross 
section, includin g en ergies outside of the region accessible experimentally at present, was 
taken from Ref. 12811 . 



Proton - proton collisions 

Rapidity distribution 

Fig. 121 shows the calculated pseudorapidity distributions for the proton - proton (antipro- 
ton) collisions. The agreement with the experimental data is quite good despite the fact that 
pp collisions present special difficulties for our approach since the value of the saturation 
momentum is rather small and non-perturbative corrections could be essential. We would 
like to point out however that the value of the saturation momentum for the proton reaches 
~ 1 GeV at the LHC energy. Our experience with RHIC data suggests that at such value 
of the saturation momentum our approach could apply with a reasonable accuracy. We thus 
expect that the very first pp data from the LHC can provide an important test of the CGC 
ideas. In Fig. 01 we plot the value of dN/drj at = as a function of energy (in this plot we 
included also the available data at lower energies). The agreement with the experiment is 
seen to be quite good. 



Total multiplicity 

Integrating Eq. (jHlll over t] in the entire region of ?7 = — In ly -i- +\nW we can calculate 
the total multiplicity in proton-proton ( antiproton) collisions. In Fig. |3 we present our 
calculation together with the experimental data taken from Refs. [4^; a good agreement 
with the data is seen in a wide range of energies. We would like to remind however that our 
predictions for the LHC energies could be as much as 1.5 times larger due to uncertainties 
in the energy behavior of the saturation scale discussed above. 
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FIG. 2: Rapidity dependence dN/dij of charged hadron multiplicities in proton - proton (antipro- 



ton) col 



Ref. 



I42 



isions as a function of the pseudorapidity at different energies. The data are taken from 



Nucleus-nucleus collisions. 



Rapidity distribution dN/dy 

Our prediction for lead-lead collision at the LHC energy is plotted in Fig. El The two 
sets of curves ( solid and dotted) describe the cases of fixed and running QCD coupling 
respectively. We consider the two predictions as the natural bounds for our predictions, and 
expect the data to be in between of these two curves. However, we would like to mention 
again that our predictions have systematic errors of about 12 15% due to uncertainties in 
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FIG. 3: Energy dependence of charged hadron multiplicity dN/drj at = in proton - proton 
(antiproton) collisions and of charged hadron multiplicities per participant pair (2 / Npart) dN / drj 
at r/ = for central nucleus- nucleus collisions. The vertical dotted lines mark the LHC energies 
for nucleus-nucleus collisions {W = 5500 GeF) and for proton-proton collisions {W = 14000 GeF). 



collisions. The experimental data are from Ref. |42 
the energy dependence of the saturation scale. 



43|. 



Centrality dependence: {2/Npart) {dNch/drj) 



Fig. El shows our predictions for the Npart dependence of the {2/Npart) [dNch/drf). This 
observable provides the most sensitive test of the value of the saturation scale and its de- 
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FIG. 4: Energy dependence of total multiplicity in proton - proton (antiproton) collisions. The ver- 
tical dotted line marks the LHC energies for proton-proton collisions (W = 14000 GeF). collisions. 
The experimental data are taken from Ref. 4^ . 



pendence on the density of the participants. 

Fig. El shows the energy dependence of {2/Nch) {dN/drj) at r/ = 0. One can see that we 
are able to describe the current experimental data. Note that if we neglect the difference 
between rapidity and pseudorapidity, {2/Nch) {dN/dif) at r/ = is given by a very simple 



formula 



id: 



dN, 



ch 



part 



dy 



\y=0 



0.87 (^) \n{Ql{AW,y = 0)/AlcD) 



0.87 



W 
130 



0.252 



(3.93 + 0.252 ln(iy/130) 



(34) 
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Lead - Lead Collision 



: dN,/dr| W = 5500GeV 




-8 -6 -4 -2 2 4 6 8 



FIG. 5: Rapidity dependence dN/drj lead- lead collisions at the LHC energy at different centrality 
cuts. The solid lines corresponds to the prediction using Eq. ()16|) for the energy dependence of the 
saturation scale while the dotted lines show the predictions for Eq. 1)26(1 for running QCD coupling. 
The shadowed area shows the prediction for the minimal bias event. 

This formula is in good agreement with the existing experimental data. 

Proton - nucleus collisions 

Fig.[7|shows our prediction for the proton-nucleus collisions at W=5500 GeV. In section 2 
we described the procedure of computing the number and density of participants in this case; 
to evaluate the relevant value of the saturation momentum, we take account of the energy 
dependence to extrapolate from RHIC to the LHC energy. For example, the density of par- 
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Centrality dependence for Pb-Pb Collision 




FIG. 6: Npart dependence of {2/Nch) dN/drj for lead-lead collisions at the LHC energy at different 
rapidity cuts. The solid lines correspond to the prediction using Eq. (jlbj) for the energy dependence 
of the saturation scale while the dotted lines show the predictions for Eq. ()26|) for running QCD 
coupling. The shadowed areas show the spread of our predictions. 

ticipants ppart ~ 1-84 fm"^ corresponds to the saturation scale of ~ 2 GeV^(5500/200)°-^^^ 
^4.6 GeV^ 

CONCLUSIONS 

In this paper we have provided a complete set of predictions for the multiplicity distribu- 
tions at the LHC basing on the CGC approach. In our approach, parton saturation results 
in a relatively weak, compared to most other approaches, dependence of the multiplicity on 
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FIG. 7: Rapidity dependence of dN/dr] proton-lead collisions at the LHC energy at different 
centrality cuts. The dotted line corresponds to the minimal bias event. 

energy. As one can see from Fig. |H1 we expect rather small number of produced hadrons in 
comparison with the alternative approaches. What is the uncertainty in our predictions? 
We would like to recall the estimates for the uncertainty in our calculations at the LHC 
energies given above: 12 -^ 15% for nucleus-nucleus and hadron-nucleus collisions, and a 
large value of 40 -t- 50% for the proton-proton collisions. These uncertainties arise from the 
poor theoretical knowledge of the energy dependence of the saturation scale, and in the case 
of pp collisions also from the uncertainties in the application of the geometrical picture. 

We hope that our estimates will be useful for the interpretation of the first results from 
LHC experiments. As illustrated in Fig. |Hl a measurement of multiplicity at the LHC will 
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FIG. 8: Comparison of our predictions for charged hadron multiplicities in central (b < 3 fm) 
Pb — Pb collisions with the results from other approaches, as given in Ref.j^l 

provide a very important test of the CGC approach. 
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